Two-component gap solitons with linear interconversion 
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We consider one-dimensional solitons in a binary Bose-Einstein condensate with linear coupling 
between the components, trapped in an optical-lattice potential. The inter-species and intra-species 
interactions may be both repulsive or attractive. Main effects considered here are spontaneous 
breaking of the symmetry between components in symmetric and antisymmetric solitons, and spatial 
splitting between the components. These effects are studied by means of a variational approximation 
and numerical simulations. 
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Introduction: The existence of quasi-one-dimensional 
(ID) solitons in Bosc-Einstcin condensates (BECs) has 
been demonstrated in well-known works which used 
the Feshbach-resonance (FR) technique to switch the re- 
pulsion between atoms into attraction. It was also pre- 
dicted 01, and then demonstrated in the experiment 3, 
that an optical-lattice (OL) potential may support gap 
solitons (GSs) in finite bandgaps of the OL-induced spec- 
trum. 

Experiments are possible too in a binary BEG, created 
as a mixture of two hyperfine states of the same atom 
[1]. In that case, the inter-species interaction may also 
be controlled by FR Q . It was proposed to use this set- 
ting for the creation of symbiotic solitons @ , in which the 
FR-induced inter-species attraction overcomes the intra- 
species repulsion. Two-component symbiotic GSs were 
predicted too Q. In particular, the attraction between 
two self-repulsive species may lead, counter-intuitively, 
to spatial splitting between GSs formed in each species 
Q , which is explained by a negative effective mass of the 
GS. A general two-component model including an OL 
potential acting on both species was considered in Ref. 
[9|, with inter-species repulsion and intra-species repul- 
sion or attraction. Several types of GSs were reported 
there: symmetric and asymmetric ones, unsplit or split 
complexes (in which, respectively, centers of the two com- 
ponents coincide or are separated) , and solitons classified 
as intra-gap or inter-gap ones, with chemical potentials 
of the two components belonging to the same or different 
bandgaps. 

Our objective here is to extend the consideration to a 
different physical situation, when the interaction between 
the two components trapped in the OL includes linear 
coupling, i.e., interconversion between the species. The 
linear interconversion between hyperfine atomic states 
can be induced by a resonant electromagnetic field p^ . 
leading to the prediction of various coupled-mode dynam- 
ical effects, such as Josephson oscillations 
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On the contrary to the previously studied symbiotic 
solitons 0, 0j B 3 , chemical potentials of the two com- 
ponents of stationary patterns in the linearly-coupled bi- 
nary BEG must be equal. This condition, obviously, fa- 
vors symmetric states, i.e., ones with identical compo- 
nents. On the other hand, nonlinear repulsion will push 



them aside. Thus, the competition between the linear 
coupling and nonlinear inter-component repulsion gives 
rise to a shift of the miscibility-immiscibility transition 
in a binary BEG, or in a superfluid Fermi gas trapped 
in a parabolic potential [l2|; however, manifestations of 
such competition in terms of solitons were not considered 
before, and this is one of objectives of the present work. 
On the other hand, the interplay of the intra-component 
self-repulsion or self-attraction and linear interconversion 
gives rise to spontaneous symmetry breaking (SSB) in ID 
[ill and 2D ^ binary solitons. The SSB manifests it- 
self in the destabilization of symmetric or antisymmetric 
solitons (in the case of the self-attraction or repulsion, 
respectively), and the emergence of stable asymmetric 
states, with different numbers of atoms in the compo- 
nents. The SSB in two-component GSs was also pre- 
dicted in optics, in terms of dual-core fiber Bragg grat- 
ings [H]. 

In this work, we use the variational approximation 
(VA) and numerical simulations to demonstrate that, in 
the presence of the linear interconversion of components, 
both the SSB and spatial-splitting transition in binary- 
BEG solitons trapped in the OL potential take a simple 
form, which can be represented by universal diagrams in 
the parameter space. 

The model: We consider a binary BEG loaded into 
"cigar-shaped" trap, combined with the OL potential 
acting in the longitudinal direction. The system of cou- 
pled Gross-Pitaevskii equations for the two wave func- 
tions, 01 and 02, can be written as [1, HBQ^ 

= -(1/2) (0,)..+ {g\q^j\^ + gi2\^3-j\^)cl)j 
-Vq cos {2x) - K(j)3^ J, j = l,2. (1) 

The scaled coordinate, time, OL strength, linear- 
interconvcrsion rate, and nonlinearity coefficients are 
related to their counterparts measured in physical 
units: x = {'K/L)xp\^, t = (tt/L)^ (/i/m) tph, Vq = 
(L/ttS)^ m(Vo)ph7 and k = (L/tt)^ (m//i) Kph, and 
{5,512} = {Lmu!±^H /i:h) {a, 012}, where L, m, w^, and 
M are the OL period, atomic mass, transverse-trapping 
frequency, and total number of atoms in both species. 
We define Vq and k to be positive, while a, ai2 > and 
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< correspond to the nonlinear repulsion and attraction, 
respectively. With t replaced by propagation distance z 
and gi2 = 0, Eqs. H]) may also be interpreted as a model 
for the spatial evolution of optical signals in a planar 
dual-core wave gui de, equipped with a transverse grating 
of strength Vb [ll. 

Stationary solutions to Eqs. ^ are looked for in the 
usual form, (j)i^2{x,t) = exp (— i^i) ui_2(a;), with a com- 
mon chemical potential, /i. Real functions ^1,2(2;) obey 
the following equations and normalization: 

^Uj+u- /2~ [gu^ + gi2ul_j) Uj+Vq cos (2.t) uj+kUj = 0, 

(2) 

A^i + A^2 = 2, Nj = u'j{x)dx, which can be derived 
from the Lagrangian, 



E 

3 = 1.2 



^lu^j - {u'jY/2 + 1/0 cos(2a;)M 



-K/2 



gi2u\u\ + 2nuiU2 ^dx — 2fj, 



(3) 



Solitons may exist if fi falls into the semi-infinite gap 
(SIG) or finite bandgaps of the linear spectrum of system 
^ , which was found in Ref . [l^ . 

The variational approximation. To predict solitons 
with a compact unsplit profile and, generally, different 
numbers of atoms in the components, A^i ^ N2, we adopt 
the Gaussian ansatz p^ . 

^unspi)^^^ = a^7r-i/4y]^g--V(2-=)^ j ^ 1^2, (4) 

where, for A^i — N2, a = +1 and —1 correspond 
to symmetric and antisymmetric solitons, respectively: 
^(unspi) _ _|_y(unspi)^ p^gg parameters in ansatz Q are 
norms A^i,2 of the components and their common width 
w (in the presence of the linear coupling, one may assume 
equal widths of the components, even if their amplitudes 
are different [H, [3). The substitution of ansatz (g]) in 
the Lagrangian yields 



L = ^i{N-2)-N/ (4u;2) + Voe' 

+ 2aKv - [g {N^ - v^) + gi2v'^] / {2y/2^w) , (5) 



where N = Ni+ N2, u = ^/WiN^. In this notation, the 
asymmetry parameter of the soliton is 



= (iVi - N2) I (A^i + N2) = v/l-i^V2. 



(6) 



The first variational equation, dL/dfx = 0, recovers 
the normalization adopted above, N ~ 2. The other 
equations, dL/dw = dL/dv = and dL/dN = 0, yield 
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2^/TrKw ^ a {gi2 - g) \/\~ t^, 



w = 0, 
(7) 
(8) 



and \i = (4w2) - Vbe"™ + y/2jTTgw-'^. Note that Eq. 
([8]) has no solutions for a (512 — g) < 0, which means that 
stationary asymmetric states with identical signs of the 
two wave functions (a = 4- 1 ) do not exist if the repulsion 
between the species is weaker than the intrinsic repul- 
sion in each of them, or, alternatively, if the inter-species 
attraction is stronger than its intra-species counterpart. 
Just the opposite is true for states with different num- 
bers of atoms and opposite signs of the wave functions, 
a = —1. The equation for width wq of symmetric and 
antisymmetric states is obtained from Eq. ([7]) by setting 
e = in it. 



WqWq exp I 



{g + g,2)/V2^, (9) 



while Eq. ([8]) should be ignored in this case, as its deriva- 
tion implied a deviation from the symmetry. 

An issue of major interest is to predict a critical value 
of the linear-coupling constant, k ~ Kbif, at which the 
SSB bifurcation occurs, i.e., solutions to Eqs. (O and 
(l8|) with infinitesimal e split off from the solutions with 
e = 0. In this case, Eq. ([7]) again reduces to but Eq. 
([8]) should not be omitted, yielding 



Kbif = cr {gi2 - g) / (2^/t:wq) . 



(10) 



The VA may also be applied predict the splitting border 
in the case of the repulsive nonlinearity. To approximate 
the onset of the splitting, we use the ansatz introduced 
in Ref. @ (in the absence of the linear coupling) , 



,(spi) 



(x) 
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l±bx- ]^v? {b^ ~ x^) 
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where the separation between the centers of the compo- 
nents is Aa; « 2b'w^ for small b. The substitution of this 
ansatz in Lagrangian ([3]) yields 
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L = -— + 2Vbe" 
-262 / 2Vow^e-'^^ 
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and the additional variational equation, dL/d{b^) = 0, 
predicts the splitting condition: 



512 



= 2\/27rVoi 



-|- a\/2TTKWo. 



(14) 



Numerical results. Numerical solutions of Eq. ([T]) were 
obtained by means of the real-time propagation using the 
split-step Crank-Nicolson algorithm Q , with spatial and 
temporal steps 0.025 and 0.0002, respectively. The sim- 
ulations were run until the solution would settle down 
into a stationary localized state. This method of obtain- 
ing stationary solutions guarantees their stability. 

The most fundamental information about the transi- 
tion between different types of two-component solitons is 
provided by the SSB and splitting borders in the plane 
of interaction coefficients {912, g), for fixed coupling con- 
stant K. A generic example of such a diagram is presented 
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FIG. 1: (Color online) Four curves on the left represent 
lines of the spontaneous symmetry breaking of symmet- 
ric ((T — +1) and antisymmetric (a = —1) solitons, for 
Vo = 5. Three curves on the right designate the onset 
of the spatial splitting between the components of symmet- 
ric and antisymmetric solitons. The results are obtained 
from the variational approximation (VA) and numerical sim- 
ulations (num). Labels "overlapping/bifurcated" and "un- 
split/split" indicate, respectively, areas of the existence of 
stable symmetric/asymmetric solitons, and of those with co- 
inciding/separated centers of the two components. 



in Fig. [l]for Vq = 5 and k = 0.5 (in physical units, the 
latter corresponds to interconversion time 50 /iS, for 
atoms of ^Li trapped in the OL with period L ^ 1 fim; 
the range of values of |5,5i2| displayed in Fig. [T] corre- 
sponds to the solitons built of up to ~ 10^ atoms). The 
variational results displayed in the diagram are obtained 
from Eqs. PT7)) and (|14[) . using a numerical solution of 
Eq. ^ for wq. In the range of 0.1 < k < 1, and for other 
relevant values of Vq, the diagram keeps essentially the 
same form as in Fig. [T] 

The diagram covers both positive (repulsive) and neg- 
ative (attractive) values of gi2 and g, the SSB line in 
quadrants (712,5 > and gi2,g < going, respectively, 
through families of GSs and regular solitons (the latter 
ones belong to the SIG). It is worthy to stress that the 
broken-symmetry areas are located on opposite sides of 
the SSB lines for a = +1 and —1, i.e., for the the sym- 
metric and antisymmetric solitons. The VA version of 
the splitting line for cr = -1-1 is not included in Fig. [U as 
in this case SSB happens prior to the onset of the split- 
ting, while Eq. (fTi|) was derived from the VA assuming 
the unbroken symmetry. Note also that the SSB and 
splitting lines do not intersect, and there is no overlap 
between the stability areas of different types of the two- 
component solitons, i.e., the model does not give rise to 
a bistability. 

The splitting lines displayed in Fig. [T] are similar to 
those reported in Ref. 0] for the model with k = (nev- 
ertheless, the difference is that k 7^ makes complete sep- 
aration of the two components of the soliton impossible, 
unlike the situation with k = 0). On the other hand, all 
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FIG. 2: (Color online) Examples of stable symmet- 
ric/asymmetric and unsplit/split solitons, as obtained from 
the numerical solution for different combinations of the at- 
tractive and repulsive interactions. The thin dashed curve 
depicts the OL potential, with Vo = 5. 



the results concerning the SSB lines have no previously 
published counterparts, as the symmetry breaking in reg- 
ular and/or gap- mode solitons was not studied before in 
models featuring the competition between the linear cou- 
pling and nonlinear interaction between the species. 

Figure[2]displays typical examples of SSB in symmetric 
and antisymmetric solitons [panels (a), (b) and (c), (d), 
respectively], and of the spatial splitting combined with 
SSB [(e), (f)]. To stress the role of the linear coupling 
in inducing the transition to the asymmetric shapes, the 
panels also include the case of k = 0. For all unsplit 
solitons, the shapes predicted by the VA, both symmet- 
ric and asymmetric ones, are virtually identical to the 
numerically found shapes shown in panels (a)-(d). 

The case when the intra-species nonlinearity is 
switched off {g = 0), and the solitons exist only due 
to the inter-species interactions, is of particular inter- 
est. The evolution of families of stable asymmetric GSs 
with the increase of [512 1, as obtained from the numeri- 
cal solution in this case (both for the attraction, 312 < 0, 
and repulsion, gi2 > 0, when the solitons can be found, 
severally, only with (7 = — lor(T = +l),is shown in Fig. 
[3l for different values of linear coupling k. In particular, 
the initial abrupt increase of e with the growth of 1712 > 
is the manifestation of the SSB in the unsplit GS. After 
achieving a maximum of e very close to 1 , the asymmetric 
GS undergoes the splitting transition, which eventually 
leads to the gradual decrease of the effective asymmetry 
between the spatially separating components. In the case 
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FIG. 3: (Color online) The asymmetry of the solitons, e, de- 
fined as per Eq. versus the nonlinear-coupling coefiicient, 
512, as found from the numerical solution for Vo = 5, g = 0, 
and cr = 1 (four lines on the right side) or cr = —1 (four lines 
on the left side). 
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FIG. 4: (Color online) The time evolution of the symmetric 
gap soliton found at cr = 1, <; = 4, 512 = 2, k = 0.2, Vo = 5, 
after as g suddenly dropped from 4 to 1. Here, (fii = (f>2 = 

4>{x,t). 



of 512 < 0, the solitons (in this case, they belong to the 
SIG), also attain a maximum of e very close to 1. Fur- 
ther increase of |(7i2l leads to a transition to solitons of 
the symbiotic type [6| with reduced asymmetry. These 
results, as well as those summarized in Fig. [U do not 
have previously published counterparts cither. 

As mentioned above, the numerical procedure adopted 
in this work generates only stable soliton solutions. How- 
ever, unlike the ordinary SIG-based solitons, GSs do not 
represent the BEG ground state 0, i-S-, they are, 
strictly speaking, mctastablc objects. For this reason, 
it is relevant to test their stability against strong pertur- 
bations. Simulations demonstrate that the GSs are, in 
fact, very robust objects. For instance, sudden drop of 
the self-repulsion coefficient from 5 = 4 to g = 1, which 
can be easily implemented by means of FR, does not de- 
stroy the GS, see Fig. U) Similarly, sudden application 
of a kick to the soliton (not shown here) gives rise to its 
oscillations around a local minimum of the OL potential, 
but does not destroy it either. 

Conclusion. In this work, we have extended the model 
of binary BEG trapped in the OL potential Q by includ- 
ing a the linear interconversion between the two species. 
Using the VA and numerical simulation, we have identi- 
fied two internal transitions in two-component solitons, 
both regular and gap- mode ones: spontaneous symmetry 
breaking, and spatial splitting between the components. 

This work was supported, in a part, by FAPESP and 
CNPq (Brazil). 
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